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ABSTRACT 

We have made a detailed spectral analysis of eleven Wide Angle Search for Planets 
(WASP) planet host stars using high signal-to-noise (S/N) HARPS spectra. Our line 
list was carefully selected from the spectra of the Sun and Procyon, and we made a 
critical evaluation of the atomic data. The spectral lines were measured using equiva- 
lent widths. The procedures were tested on the Sun and Procyon prior to be being used 
on the WASP stars. The effective temperature (Teff), surface gravity (logg), microtur- 
bulent velocity (wmic) and metallicity were determined for all the stars. We show that 
abundances derived from high S/N spectra are likely to be higher than those obtained 
from low S/N spectra, as noise can cause the equivalent width to be underestimated. 
We also show that there is a limit to the accuracy of stellar parameters that can be 
achieved, despite using high S/N spectra, and the average uncertainty in Tcs, logg, 
Vmic and metallicity is 83 K, 0.11 dex, 0.11 kms~^ and 0.10 dex respectively. 

Key words: stars: abundances, stars: fundamental parameters 



1 INTRODUCTION 

Since the first exoplane t was discovered orbiting 51 Pegasi 
jMavor fc Quelodl 19951 ). hundreds more have been detected 
around other stars. Due to biases inherent in the detection 
techniques, the majority of the planets discovered to date are 
giant planets with short orbital periods. How ever, with th e 
advent of the Kepler space mission (Borucki et al.ll20lj ). 
planets with properties more akin to Earth are now be- 
ing discovered. Accurate planetary parameters are essential, 
whether one is trying to study hot Jupiters or find another 
planet like Earth. In the case of transiting planets, informa- 
tion such as the orbital period can be extracted directly from 
the light curve, however the mass and radius of the planet 
are coupled with the mass and radius of the host star. The 
transit of an exoplanet across its host star will only yield the 



* The analysis is based on data obtained through observing pro- 
grams 072.C-0488, 082.C-0040, 082.C-0608, 283.C-5017, 084.C- 
0185 & 087.C-0649, from ESO's HARPS spectrograph, mounted 
at the 3.6 m at La Silla, Chile, 
f E-mail: a.doyle@keele.ac.uk 



ratio of the planet to star radius and the mass of a planet, 
acquired from combining transiting and radial velocity data, 
cann ot be determ ined independently from the mass of the 
star llWinnlbOld ). 

It is thus imperative that the mass and radius of the 
host star are known precisely, but obtaining direct measure- 
ments is only possible for a limited number of stars. For in- 
stance, a fundamental value of mas s cannot be obtained un- 
less t he star is in a binary system (|Torres. Winn Sz HolmanI 
|2008| ). and determining the radius of a star requires knowl- 
edge of t he angular diamete r, which in turn needs a known 
distance l| North et al.l |2007| ). Most current distance mea- 
surements were acquired using the parallax technique with 
ESA's Hipparcos satellite, and these will be complemented 
in the near futu re by data from ESA's Gaia mission, due to 
launch in 2013 jCare. Sarri fc Schmidti,20ia '). 

In the absence of direct measurements, stellar spec- 
troscopy can be utilised to determine the effective tem- 
perature (Tcff), the surface gravity (logg), and metallicity. 
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[Fe/HB of the planet host star. These parameters are then 
used to infer the stellar mass and radius, based on a calibra- 
tions such as [Torres. Anderson k, Gimenej ||2010| ) , or a grid 



Table 1. Details of the HARPS spectra used in the present work 



of stellar models such as Girardi et al.l ( 2000l ). 

Knowledge of the chemical composition of the host 
stars has additional value, as planets are more likely to 
be found around high metallicity stars jGonzalej 19981 : 
ISantos. Israehan fc Mayorll2004l : [Fischer fc Valentill2005l ). It 
has also been suggested that stars with planets can have dif- 
ferent chemical compo sitions to stars without known planets 
l|Melendez et al.ll2009l ). 

In order to extract the maximum amount of informa- 
tion from the statistics of planet frequency and stellar abun- 
dances, these abundances need to have been measured with 
the highest level of precision possible. However, determin- 
ing abundance is by no means an exact process. They are 
heavily influenced by the reliability of atomic data and the 
quality of the spectra. 

Details about the electron transitions that cause spec- 
tral lines, along with data on how these lines can be broad- 
ened by atomic processes, ideally need to be obtained via 
laboratory measurements. However, the sheer multitude of 
spectral lines makes this a lengthy and difficult task, so that 
often atomic data is inferred from indirect methods. The re- 
sult of this is that even spectra with a high signal-to-noise 
ratio (S /N) there is a limit as to how precisely one can mea- 
sure elemental abundances in stars. 

Spectral analysis of low S/N spectra is no easy task, as 
the excessive noise in the continuum of these spectra makes 
normalisation difficult, causes the wings of strong lines to be 
underestimated, and depletes the number of weak lines avail- 
able for measurement. The initial analyses for the WASP 
planet host stars for the discovery papers are usually per- 
formed using coadded spectra from the CORALIE spectro- 
graph, which has a resolution betwee n 55 000 and 60 000 
(Oueloz et al]|200ll : lw"ilson et al.ll2008l '). These spectra typi- 
cally have S/N varying from 50:1 to 100:1. The analyses per- 
formed here use coadded spectra obtained with the HARPS 
spect rograph which has a resolution of 115 000 (|Mavor et al.l 
l2003h . and produces spectra with high S/N, as shown in Ta- 
ble [1] Thus these spectra are more suitable for determining 
stellar parameters. The spectra were reduced using the stan- 
dard HARPS data reduction software and the observatio n 
details for these spe ctra are discussed i n iTriaud et aLll|2010^ ■ 
iQueloz et all (|2010l ) and lGiUon erall l|2009l ). 

The paper is divided as follows: section 2 outlines the 
methods used in the analyses, section 3 gives the parameters 
of the standard stars that were used to test the methods, sec- 
tion 4 details a discussion of the results, and the conclusion 
is given in section 5. 



Star 



V mag No. of spectra S/N 



^ Abundances are given in the format of log(yl) 12, where 
log(A) is the number ratio of the element with respect to hydro- 
gen, log (N^i/Nh)- When square brackets are used, e.g. [Fe/H], 
it denotes the abundance relative to the solar value, where as 
logyl(El) indicates log{A) + 12. Using curved brackets, (Fe/H), 
can also denote log(A) + 12 



WASP-2 


11.98 


32 


100 


WASP-4 


12.5 


32 


100 


WASP-5 


12.3 


34 


125 


WASP-6 


11.9 


44 


125 


WASP- 7 


9.5 


20 


235 


WASP-8 


9.79 


80 


430 


WASP- 15 


11.0 


51 


240 


WASP-16 


11.3 


77 


175 


WASP- 17 


11.6 


60 


230 


WASP- 18 


9.3 


21 


270 


WASP- 19 


12.59 


36 


125 



2 METHOD 

The spectral synthesis package used was UCLSYN ( U niver- 
sity College London SYNthesis: iSmith fc Dworet skv' 198^; 
ISmithI 1 19921 : ISmaUev. Smith fc Dworetskvl l20oiI l ATLAS9 
mode ls atmospheres without convectiv e overshooting are 
used l|Castelli. Gratton fc Kurucz|[T997l ') and local thermo- 
dynamic equilibrium (LTE) is assumed. 



2.1 Line list 

It is impractical to individually measure every line in a spec- 
trum, and so an appropriate selection of lines must be chosen 
from a high S/N sp ectrum, such as the Kurucz solar atlas 
(|Kurucz et al.|[T984l ) as was used here. Lines were also se- 
lected from the HARPS Procyon (HD 61421) spectrum, so 
that these could be used for hotter stars in addition to the 
main line list. The line list was constructed by prioritising 
unblended lines where possible, and these comprise around 
one third of the list. Some lines had blending components 
that could be resolved via spectral synthesis, and these were 
included to increase the number of lines. Unresolved blends 
were included when there was a paucity of lines for a partic- 
ular element. In this case, a restriction was imposed to en- 
sure that the blending component comprised no more than 
five per cent of the overall equivalent width (EW), other- 
wise the line was deemed to be too severely blended to be 
of use. Lines which were blended with the same element in 
a different ionisation state were rejected, even if the weaker 
component was less than five per cent, as lines blended with 
other ionisation states could bias our estimate of the Tctt- 
Unresolved blends were rejected for Fe, as there are still a 
sufficient number of Fe lines remaining after these blends are 
excluded. It is important to include as many Fe ll lines as 
possible, as these are essential when determining logg via 
the ionisation balance method. 

A multitude of lines will not necessarily decrease the 
abundance error, espe cially if the add itional lines are of poor 
quality. According to iKurucj (|200!j h , abundance errors are 
likely to increase if there is a wide range of line strengths in 
a line list. An accurate abundance can theoretically be de- 
termined from a single weak line, as weak lines (less than 30 
mA) are less affected by damping parameters and microtur- 
bulence (umic). As such, including strong lines can increase 
the abundance errors. However, it is not possible to use only 
weak lines as there are not a sufficient number of them, and 
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Table 2. A sample of the line list used in the analyses of the 



HARPS spectra. 


Element 


Wavelength 


X 


logs/ 


NIST rating 


Co 11 


4516.633 


3.459 


-2.562 




Fe II 


4546.467 


4.186 


-2.510 




Ca I 


4578.551 


2.521 


-0.697 




Ti II 


4583.409 


1.165 


-2.870 


D 


Cr II 


4588.199 


4.071 


-0.627 


D 


Cr II 


4592.049 


4.074 


-1.221 


D 


Fe 11 


4620.521 


2.828 


-3.210 


C 


Fe I 


4631.486 


4.549 


-1.594 




Fe II 


4656.981 


2.891 


-3.600 


C 


Fe II 


4720.149 


3.197 


-4.480 


D 


Mn I 


4739.087 


2.941 


-0.490 


B 



The NIST rating is given for each line where available. The 
complete table is available as supplementary content in the 
online version of this paper. 

thus strong lines were also included in the list. The EWs for 
the Fe I lines in the HARPS solar spectrum range from 5.6 to 
133.1 mA (although the strongest lines were often rejected 
at a later stage) and the lower level excitation potential (x) 
ranges from 0.052 to 5.033 eV. 

Atomic data w as mainly taken from the VALD database 
l|Kupka et al.l 19991). however dam ping parameters were ob- 
tained from iKurucz fc Beil l|l995l ) when no data was avail- 
able in VALD. In addition, the default oscillator strength 
{log gf) value was not always used for certain elements 
(e.g. Cr II, Ti ll), as the value given resulted in abun- 
dances that were highly inconsistent with the assumed solar 
abundances. For instance, using the default VALD values 
yiel ds a log^(Cr l l ) of 5 .80, which is 0.16 dex higher than 
the lAsplund et al] (|2009| ) abundance. However, using alter- 
native logg/ values from VALD results in a logj4(Cr ll) 
of 5.66. In these cases the different source of log gf from 
within VALD was used, so that a more suitable solar abun- 
dance could be obtained. The log gf values were also sup- 
plern e nted from other sources wh ere possible ijFuhr fc Wiesd 
I2OO6I : iMelendez fc BarbuvllioogI '). The choice of log 5/ can 
greatly influence the abundance obtained, resulting in a 
wide range of abundances obtained for a particular ele- 
ment, even in the Sun. For example, the solar Fe abun- 
dance varies between 7.4 1 and 7.56 throughout the literature 
(jMashonkina et aLll201ll), and the solar Mn ranges between 
5.23 and 5.46 l|Bergemann fc Gehrenll2007l ). 

All lines were cross referenced with the NIST database 
to check the reliability of log gf values. Any lines with an 
"E" ratin g, i.e. with an uncertainty of greater than ± 50 
per cent (jFuhr fc Wiesd [20061 ) . were rejected. However, it 
should be noted that not all lines in the list were present 
in the NIST database, so some lines that were included in 
our line list could still have large uncertainties in the log gf 
values. The line list is given in Table [J] showing the NIST 
rating where available. 



2.2 Balmer lines as temperature indicators 

The Balmer lines can be u sed a s an i ndication of T^s for 
stars cooler than 8000 K (jCravl 1200^ . The core to wing 
transition is the part of the line which is most sensitive to 
temperature changes l|Mashonkinall2010l ) and the outer wing 



broad ening of Hg is barely affected by other stellar param- 
eters IStempels et al.ll2007^ ■ However, Balmer lines are not 
useful for determining temperature in the HARPS spectra, 
as the over lapping echelle order s are merged during pipeline 
reduction (jBruntt et al.ll2010al ). This affects Ha more than 
H/s, and makes normalisation extremely difflcult. As such, 
the Balmer lines were only used to gain an initial estimate 
for r,ff. 



2.3 Surface gravity from pressure-sensitive lines 

Strong, pressure-broadened lines in the spectra were fitted 
in order to obtain a value for log g. The Na I D lines were 
the main lines used, however it can be difficult to acquire 
the logg from these lines if interstellar Na is present. The 
abundance used can also influence the loggr measurements. 
For example, there are only two suitable unblended Na lines 
in the spectra, other than the Na I D lines, which can be used 
for abundance determination. The lack of lines could give 
an erroneous abundance, which is then incorporated into 
the log g determination. There are more Ca lines which can 
be used to determine abundance, which led to the Ca line 
at 6439 A also being used as a gravity diagnostic. However 
this line is not as useful in hotter stars, as the sensitivity to 
changes in logg is reduced. The Mg I b lines were not used, 
as the only line that gives a reasonable value for log g is the 
5167 A line, and this is severely blended. 



2.4 Line Broadening 

Microturbulence and macroturbulence (umac) are ad hoc 
broadening parameters that are needed in order fit the 
line profiles to the observed spectra. It was found that the 
abundance derived from strong lines differed from those 
calculated using weak lines, and thus i ntroducing fmic 
woul d bring the abundances into agreement (jStruve fc ElvevI 
I193 4I). Microturbulence represents a small scale Gaussian 
velocity field, that is in added to the thermal Doppler 
broadening. It results in strong lines being spread over 
a larger wavelength interval and beco ming saturated a t 
larger EWs than if no Vmic was included (|Landstreetll2006l ) . 
Macroturbulence only broadens the lines and does not 
change their line strengths. The size of the microturbu- 
lent cell is defined as being less than the mean free path 
of the photon, whereas Wmac represents velocities that oc- 
cur where the cell is larger than the unit optical depth 
(|Howarth|[2003l ). This division between micro and macro- 
turbulence is crude, as in reality the motion would occur 
over a range of scales (iMucciarelli I2OIII : lAsplund et al.ll200ol : 



iNordlund. Stein fc Aspltmd 



2009f) . It has been shown that 
3D models have no need for Vmic and t^mac, as c onvective 
energy transport can be included in the models (jAsplundl 
l2005h . 

Microturbulence was determined from reducing the 
trend in EW with abundance. Weak lines should not be af- 
fected by Wmic, therefore requiring that the abundances of 
strong lines equal those o f the weak lines should result in a 
value for Umic- While the iMagainI l|l984h method of deter- 
mining Vmic using calculated EWs agrees with the observed 
EWs (within the errors) for spectra with S/N greater than 
70 (.MucciarelU ,2011i ') , there is still a small differ ence be- 
tween the two methods for some stars. As such, the iMagainI 
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l|l984l method was preferred for calculating the Umic- The 
scatter in abundances makes it difficult to get get a preci- 
sion in Umic any better than 0.08 kms~^. Also, obtaining 
the appropriate Vmic value depends strongly on knowing the 
correct value of Teg and loggr. If one of these two parame- 
ters is incorrect then the iimic will be skewed. For instance, a 
change in Tcff or logg by 100 K or 0.1 dex respectively will 
result in a change in Umic of around 0.1 kms~^. 

Macrotu rb ulence was calculated based on the 
iBruntt et all (|2010al ') calibration for stars with Tch 
below 6500 K, and this was extrapolated for stars greater 
than 6500 K. Rotational velocity (usini) was then fit for 
a selection of unblended Fe I lines. The radiative damping 
constant. Van der Waals damping constant, and the Stark 
broadening factor were input via the line list. The instru- 
mental broadening, determined from the telluric lines at 
around 6300 A, was also included. 



2.5 Parameters from Fe lines 

The T^s , log .9 and fmic of a star can be determined using 
Fe lines jGonzalez fc Vanture'l998':'Brunt t. De Cat fc Aertd 
[sin 



Microturbulence was determined by requiring a null de- 
pendence between Fe abundance and EW (see Sect. 12. 4( ). 



I2OO8I : [San tos. Israclian & Mayor 2000; S ousa et all l2006h 



Theoretically, this can be done with any element, but only 
FGK stars have a suitable number of Fe lines present to 
perform a precise analysis. 

Fe I abundances will increase with increasing T^s- This 
temperature sensitivity is greater for low x lines, and is al- 
most negligible for high x lines. Thus requiring that there 
is no trend between x and abundance should yield the Tcff 
of the star. The same principle can be applied to Fe 11 lines, 
where in this case it is the high x lines that are sensitive 
to changes, however there are not usually a sufficient 
number of Fe 11 lines present. The Fe 11 abundance will in- 
crease with increasing log g, where as the Fe I abundance is 
insensitive to log;? variations. Thus, requiring that the Fe I 
and Fe 11 abundances agree should result in a value for log g. 
These methods for determining the Tcfi and \ogg from Fe 
lines a re known as excitation and ionisatio n balance respec- 
tively (|Takeda. Ohkubo fc Sadakanell2002l ). 

The mean of the logg determined from the ionisation 
balance and from the pressure-broadened lines was used 
as the overall value. The logg from the ionisation balance 
agreed with the \ogg from the pressure-broadened lines to 
within 0.3 dex. This difference is most likely due to the er- 
rors in abundance, discussed in Sect. 12.31 or due to T^s er - 
rors influencing the ionisation balance. IBruntt et al.l (|2010bl ) 
also recommend the use of the Ca I 6122 and Ca I 6162 A 
lines, however we found that these lines often gave a logg 
value that differed greatly from the logg determined from 
the other pressure-broadened lines and the ionisation bal- 
ance, therefore these lines were not included in the overall 
logg estimation. 

Any obvious outliers in the plot of Fe abundance against 
EW or X were removed, with a requirement that the abun- 
dance of the lines should be within 0.25 dex of the aver- 
age abundance. This restriction allowed discrepant lines to 
be removed while still retaining a sufficient amount of lines 
needed to determine the stellar parameters. In general it 
was found that the same lines had abundances that were 
too high or too low, indicating that these lines had incorrect 
atomic data. 



2.6 Line measurements and continuum 
determination 

We calculated the abundances by using UCLSYN to inter- 
actively measure the EWs of as many absorption lines as 
possible for each element. A direct numerical integration of 
the line profile was performed in order to determine the EW. 
A synthetic profile was then generated and compared to the 
observed profile shape. If the synthetic profile failed to pro- 
duce an acceptable fit to the observed spectrum, for example 
the observed spectral line could have been broader than the 
synthetic spectral line, then a least squares method fit was 
used to obtain the EW. 

In order to directly measure the EWs, a point on each 
side of the line is selected where the wing of the line reaches 
the continuum. However determining these points can be 
difficult, particularly for strong lines, and often results in an 
underestimate of the EW and thus the abundance. 

Uncertainties in continuum placement can also affect 
the measurement of EWs, and thus influence the abundance. 
The continuum was carefully normalised by eye over a small 
wavelength range for each line in order to ensure maximum 
precision. It was found that in spectra with low S/N (such 
as the CORALIE spectra), the noise makes the continuum 
placement difficult. In addition, the wings of the lines be- 
come lost in the noise, leading to an underestimate in the 
line strength. An example is shown in Fig. [T] for a line in 
the CORALIE WASP- 16 spectrum (S/N of 70), as well as 
for the same line in the HARPS spectrum which has S /N of 
175 

A high S/N spectrum doesn't necessarily eliminate all 
problems associated with determining where the line wings 
meet the continuum. Extremely weak lines that are ordinar- 
ily lost in the noise at the continuum become evident, but 
these are often unidentified lines which makes them difficult 
to synthesise, as is shown in Fig.[2]for the Kurucz solar atlas. 



3 STANDARD STARS 

Before applying the method outlined above to the HARPS 
spectra of the WASP stars, it was important to verify that 
the method achieved acceptable results on well known stars 
such as the Sun and Procyon. Procyon was chosen as a hot- 
ter comparison star than the Sun, as many WASP stars have 
temperatures higher than the solar value. The parameters 
were derived from the Kurucz solar atlas (S/N ~3000) and 
the HARPS sky spectrum (S/N ~1000) for the Sun, and 
from a HARPS spectrum of Procyon. These parameters, 
along with literature values, are displayed in Table [3] The 
parameters of the Sun are well known, however the param- 
eters of Procyon are less accurate. The mass, and thus the 
logg, can be determined with a good deal of accuracy due 
to the binary nature of Procyon. However, there is still some 
disagreement as to the TeS, mainly due to different values 
of bolometric flux and an gular diameter. The Tes values 
vary from 6516 ± 87 K jAu fdcnber g. Ludwig fc Kervellal 



vary irom boib ±0/ is. JAu idcnber g. IjUdwig Kervellal 
[20051), 6530 ± 49 K jAUende Prieto et al.l I2OO2I') and 



6591 ± 43 K jChiavassa et all I2OI2I '). ICasagrande et all 
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Figure 1. A Ni I line in a CORALIE spectrum of WASP- 16 with 
S/N of 70 is shown in the top panel. The noise in the continuum 
makes it difficult to normalise the spectrum, as well as creating 
problems with fitting the wings of the line. The same line is shown 
in the bottom panel with the HARPS spectrum (S/N 175) for 
comparison. 



also obtained a value of 6626 ± 80 K using the 
Infrared Flux Method (IRFM; iBlackwell fc Shaiiii Il977l : 
iBlackwell. Petford fc Shallislll980l ). An averaged T^s value 
is adopted for Table [S] 

In this analysis the solar logyl(Fe) was determined to 
be 7.52 ± 0.08 (for the HARPS solar spectrum). This is in 
good agreement with the va l ues of 7.50 ± 0.04 and 7.52 ± 
0.06 found bv Asplund et al.1 (|2009l ') and lCaffau et al.l (120111 ) 
respectiv ely, but higher than the valu e of 7.45 ± 0.02 deter- 
mined bv lMelendez fc Barbuvl l|2009l ). 



4 RESULTS AND DISCUSSION 

The updated T^a, logg and log A(Fe) values for the WASP 
stars are given in Table]?] along with the initial spectroscopic 
TeS, \ogg and logA(Fe) from the discovery papers. These 
parameters were obtained mainly from the CORALIE spec- 
tra, except for WASP-2 which used a SOPHIE spectrum and 
WASP-8 which also used a HARP S spec trum. The IRFM 
Tea from lMaxted. Koen fc SmallevI (|201ll ) is also included, 
and the T^b derived from the HARPS spectra is compared 
with the IRFM TcS in Fig.O A 1:1 relationship shows that 
there is a good agreement between the two different meth- 
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Figure 2. The top panel shows the same Ni I line as Fig.[T] but in 
the Kurucz solar atlas (S/N of 3000). The lower panel is a close 
up of the continuum showing that weak lines can cause issues 
even in high S/N spectra. 



ods, with the average T^b difference being 31 ± 110 K. How- 
ever, it is possible that the spectroscopic temperatures are 
higher than the IRFM values, particularly as the T^b we de- 
termined for Procyon is higher than the average literature 
value. In the case of WASP-17, strong interstellar reddening 
is present which could affect the IRFM T^b- 

Table [S] gives the Wmic, I'mac and iisini. The mass and 
radius for each star are obtained from the Tcb, logg, and 
[Fe/H], based on the [Torres et al.1 l|2010l ) calibration. How- 
ever, it should be noted that these masses and radii are only 
given as an example, and that stellar evolutionary models 
may give different results based on the same spectroscopic 
parameters. 

The derived abundances are given in Table [S] for all el- 
ements that have three or more usable spectral lines. The 
abundances obtained from the HARPS so l ar sp ectrum are 
largely consistent with the lAsplund et al.1 (|2009| ) solar val- 
ues. However, there is a discrepancy with Co, and the er- 
ror on Mn is large compared to the other elements. This 
may be due to hyperfine splitting as discussed in Section 
Overall, it was found that the logyl(Fe) values derived from 
the HARPS spectra were an average of 0.09 ± 0.05 dex 
higher than the previous analyses. This is excluding WASP- 
2, WASP-4, and WASP-5, as there are large uncertainties in 
the initial values. In addition, the initial analysis of WASP- 
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Table 3. Parameters obtained for the Sun and Procyon. 



Parameter 


Sun (Kurucz) 


Sun (HARPS) 


Sun (literature) 


Procyon (HARPS) 


Procyon (literature) 


Teff(K) 


5760 ± 50 


5775 ± 45 


5777" 


6660 ± 95 


6566 ± 65'' 


logg 


4.42 ± 0.02 


4.43 ± 0.02 


4.44° 


4.05 ± 0.06 


4.01 ± 0.03^= 


logyl(Fe) 


7.49 ± 0.06 


7.52 ± 0.08 


7.50 ± 0.04'* 


7.48 ± 0.09 


7.36 ± 0.03'= 


iJmic (kms-1) 


0.85 ± 0.08 


0.75 ± 0.15 


0.85^ 


1.70 ± 0.08 


2.2<= 



°jGr31iQQi), See text, '^ IChiavassa et al.l l l2012h . '' lAsplund et al.l | |2009|) . " lAUende Prieto et al.l ||2002| '). J' iMagainI 
||1984| 1 
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Figure 3. Comparison of T^ff from the HARPS spectra with the 
IRFM. The solid line depicts the 1:1 relationship. 



8 was also of a HARPS spectrum, so that the abundances 
agree between both examinations of the spectrum. Higher 
abundances are to be expected from the high S/N HARPS 
spectra, as the line profile wings are less likely to be under- 
estimated due to noise on the continuum, however a detailed 
comparison is still difficult until more HARPS spectra have 
been analysed. 

In addition to using abundances determined from the 
EWs measured using UCLSYN, a differential analysis with 
respect to the Sun and Procyon was also performed. Both a 
line by line differential analysis was used, as well as an al- 
ternate method where the solar log A{Pe) was fixed, and the 
loggf values were adjusted accordingly. It was found that 
the UCLSYN abundances agree with the differential abun- 
dances to within 0.04 dex, and that the scatter is reduced 
for the differential abundances. However, a differential anal- 
ysis is best performed when the stars have almost identi- 
cal p arameters to that of the reference star (iTakeda e t al.l 
|2005| ). and adjusting the solar log gf values assumes that an 
accurate solar log^(Fe) is known. As such, the abundances 
obtained from the non-differential analyses were retained. 



4.1 Errors 

Errors in Toff were calculated from the 1-a variation in 
the slope of abundance against excitation potential, and 
range be tween 70 and 105 K. These values are consis- 
tent with [Torres et ahl l|2008l ). who suggest that TcB errors 
should not fall below 50K, despite the fact that some au- 
tomated spectroscopic analyses often give errors that are 
much lower than this. To justify this decision they cite, for 



example, the difference of around 100 K between excita- 
tio n equilibrium measurements and the IRFM determ ined 
by iRami'rez fc Melended l|2004h . iMaxted et al] (|201ll ) also 
performed a comparison between spectroscopic methods and 
the IRFM, which supports that temperature errors should 
not be any lower than 50 K. 

Variations in T^s, log g and Umic can affect the elemental 
abundances. Table[7]and Table[8]list the uncertainties when 
the stellar parameters are varied by their average errors for 
the HARPS solar spectrum and Procyon respectively. Ele- 
ments such as V are more sensitive to changes in temper- 
ature than others due to a restricted range in excitation 
potential. Therefore a large error on T^s will give V a larger 
error than the other elements. For certain elements, addi- 
tional sources of uncertainty have to be considered. Abun- 
dances can be overestimated for Mn and C o, as the line 

§5files are altered due to hyperfine splitting (jSchuler et al.l 
11*). 

A change in Tcff will affect the ionisation balance, for 
example an increase in T^s o f 100 K will result in log g in- 
creasing by 0.1 dex l|Brunttl 12009. ). This was accounted for 
by varying the T^g by 1-a when using the ionisation balance 
method. The loggr determined from the pressure-broadened 
lines is also dependent on the abundance, as an increased 
abundance will cause the line to be stronger. The abun- 
dances were thus varied by 1-a when determining log g from 
fitting the Na I D and Ca I lines. 

While there is some reduction in errors of stellar pa- 
rameters from the original low S/N spectra, the errors are 
still high considering the S/N of the HARPS spectra. We 
calculated the average random errors in Toh, logg, Vmic and 
logA(Fe) to be 83 K, 0.11 dex, 0.11 kms"^ and 0.10 dex 
respectively for this sample of stars. The scatter in Fe abun- 
dances due to uncertainties in atomic data can influence the 
log A(Fe) as well as T^g, \ogg and Umic- 

Systematic errors in the EWs can be investigated by 
comparing the E Ws of the s ame star between two differ- 
ent spectrographs l|Gratton et al. 2007). Figure U shows the 
EWs of the Sun measured from the Kurucz solar atlas plot- 
ted against the EWs from the HARPS solar spectrum, along 
with a 1:1 relationship. The EWs agree between both spec- 
tra for weak lines, however the HARPS EWs seem to be 
weaker than the Kurucz solar atlas for stronger lines. There 
is no significant offset evident when comparing the results 
between the HARPS solar spectrum and the Kurucz solar 
atlas, as seen in Table [S] In addition, the strongest lines are 
usually culled prior to the final analyses. Therefore any sys- 
tematic errors present should have a negligble effect on the 
stellar parameters. 
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Table 4. Results from HARPS spectra compared with the initial analyses and the IRFM T^ff. 



7 



star HARPS T;,ff (K) Initial T^ff (K) IRFM T„ff (K) HARPS log 9 Initial log g HARPS log A(Fc) Initial log A (Pe) 



WASP-2 


5175 ± 95 


5200 


± 


200" 


5110 


± 60 


4.46 ± 0.12 


4 


3 ± 0.3" 


7.46 


± 0.10 


~7.54" 


WASP-4 


5400 ± 90 


5500 


± 


150*' 


5540 


± 55 


4.47 ± 0.11 


4 


3 ± 0.2'' 


7.42 


± 0.13 


7.54 ± 0.20'' 


WASP-5 


5690 ± 80 


5700 


± 


150= 


5770 


± 65 


4.28 ± 0.09 


4 


3 ± 0.2= 


7.63 


± 0.10 


7.54 ± 0.20= 


WASP-6 


5375 ± 65 


5450 


± 


lOO'' 






4.61 ± 0.07 


4 


6 ± 0.2^* 


7.35 


± 0.09 


7.34 ± O.IO"* 


WASP-7 


6550 ± 70 


6400 


± 


100= 


6520 ± 70 


4.32 ± 0.06 


4 


3 ± 0.2 = 


7.68 


± 0.06 


7.54 ± 0.10= 


WASP-8 


5560 ± 90 


5600 


± 


so-f 


5570 


± 85 


4.40 ± 0.09 


4 


5 ± 0.1-f 


7.70 


± 0.11 


7.71 ± 0.07^ 


WASP-15 


6405 ± 80 


6300 


± 


1009 


6210 


± 60 


4.40 ± 0.11 


4 


35 ± 0.159 


7.52 


± 0.10 


7.37 ± 0.11" 


WASP-16 


5630 ± 70 


5700 


± 


150'' 


5550 


± 60 


4.21 ± 0.11 


4 


5 ± 0.2'' 


7.59 


± 0.10 


7.55 ± O.IO*" 


WASP-17 


6700 ± 105 


6550 


± 


100* 


6500 


± 75 


4.34 ± 0.23 


4 


2 ± 0.2* 


7.40 


± 0.10 


7.29 ± 0.09' 


WASP-IS 


6400 ± 75 


6400 


± 


100^" 


6455 


± 70 


4.32 ± 0.09 


4 


4 ± 0.155" 


7.60 


± 0.08 


7.54 ± O.O9J 


WASP-19 


5460 ± 90 


5500 


± 


100*= 


5440 


± 60 


4.37 ± 0.14 


4 


5 ± 0.2*= 


7.66 


± 0.11 


7.56 ± 0.09*= 



"ICollier Cameron et al.l i2007l). ''IWilson et al.l 1200^). =lAnderBon et all i2008|l. ''iGillon et al.l )2009ll . = iHellier et al.l i2009l '). ^ iQueIoz et al.l 
)2010l ). « IWeat et al.l )2009l) . " iLister et al.l i2009l ). ' lAnderson et al.l i2010D . ^ iHellier et al.l i2009bD . iHebb et al.l )2010li 



Table 5. Stellar parameters for the WASP stars using the HARPS spectra. 



Star 



Dniic (kms 1) 



(kms vs'mi (kms ■"-) Mass (Mq) Radius (Rq) 



WASP-2 





70 


± 


0.15 


0.9 


± 


0.3 


1.9 ± 0.7 


0.87 


± 


0.07 


0.90 


± 


0.14 


WASP-4 





85 


± 


0.10 


1.4 


± 


0.3 


3.4 ± 0.3 


0.92 


± 


0.07 


0.92 


± 


0.13 


WASP-5 





75 


± 


0.10 


2.2 


± 


0.3 


3.9 ± 0.2 


1.10 


± 


0.08 


1.24 


± 


0.15 


WASP-6 





70 


± 


0.20 


1.4 


± 


0.3 


2.4 ± 0.5 


0.87 


± 


0.06 


0.77 


± 


0.07 


WASP-7 


1 


40 


± 


0.08 


5.2 


± 


0.3 


18.1 ± 0.2 


1.34 


± 


0.09 


1.32 


± 


0.11 


WASP-8 





95 


± 


0.15 


1.9 


± 


0.3 


2.7 ± 0.5 


1.04 




0.08 


1.05 




0.12 


WASP-15 


1 


15 


± 


0.08 


4.6 


± 


0.3 


4.9 ± 0.4 


1.23 


± 


0.09 


1.15 




0.16 


WASP-16 





85 


± 


0.10 


2.1 


± 


0.3 


2.5 ± 0.4 


1.09 


± 


0.09 


1.34 


± 


0.20 


WASP-17 


1 


40 


± 


0.10 


5.8 


± 


0.3 


9.8 ± 1.1 


1.29 


± 


0.12 


1.27 


± 


0.38 


WASP-18 


1 


15 


± 


0.08 


4.6 


± 


0.3 


10.9 ± 0.7 


1.28 


± 


0.09 


1.29 




0.16 


WASP-19 


1 


00 


± 


0.15 


1.6 


± 


0.3 


5.1 ± 0.3 


1.01 


± 


0.08 


1.07 


± 


0.19 



The stellar mass and radi us are estimated usin g the iTorres et al. I bold ) cahbration, and 
'J^mac is calculated from the iBruntt et al. I ll2010ai') calibration. 



.< 
E 



LU 

E 



120 



100 - 




< 
I 



40 60 80 

Kurucz solar atlas EW (mA) 

Figure 4. Comparison of EWs measured in two different solar 
specta. The solid line depicts the 1;1 relationship. 




6800 6600 6400 6200 6000 5800 5600 5400 5200 5000 

7-el, (K) 

Figure 5. Microturbulence is shown to increase with temper- 
ature. The Bruntt et al. (2010a) calibration is marked by the 
curved solid line, while the Valenti & Fischer (2005) fixed value 
of 0.85 kms~^ is represented by the horizontal dashed line. 



4.2 Microturbulence variations with T^a 

The Dmic was found t o increase with Tch, in a greem ent with 
iBruntt et all (|2010al '). [Landstreet et al] (|2009l '). and lSmallevI 
1 20041 ). It has been sugseste d that f mic should be fixe d to the 
solar value of 0.85 kms"^ (| Valenti fc Fischer] lioosl ) . How- 
ever, while this assumption is valid for solar-like stars, it is 
not appropriate for hotter stars, as is shown in Fig. [5l and 
this could skew other stellar parameters for hotter stars. 



5 SUMMARY AND CONCLUSION 

In an effort to improve the precision of stellar parameters 
of eleven WASP stars from the initial analyses, we have ob- 
tained parameters using high S/N HARPS spectra. We cre- 
ated a new line list using the VALD database, but also sup- 
plementing atomic data from other sources when necessary. 
The line list was created by selecting lines from the Kurucz 
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Table 6. Abundances for the HARPS solar spectrum and the WASP stars. 



Star 


log 


A(Ca) 


log 


A{Sc) 


log 


A(Ti) 


log 


A(V) 


log 


A(Cr) 


log A{Mn) 


log 


A{Co) 


log 


A{Ni) 


log 


A{Y) 


Sun 


6.33 


± 


0.09 


3.10 


± 


0.07 


4.95 


± 


0.08 


3.92 


± 


0.06 


5.65 


± 


0.06 


5.50 


± 


0.14 


4.84 


± 


0.11 


6.24 


± 


0.06 


2.24 


± 


0.07 


Procyon 


6.37 


± 


0.09 


3.12 


± 


0.12 


4.98 


± 


0.10 


3.93 


± 


0.03 


5.65 


± 


0.11 


5.23 


± 


0.06 


4.83 


± 


0.10 


6.20 


± 


0.08 


2.28 


± 


0.07 


WASP-2 


6.36 


± 


0.08 


3.15 


± 


0.07 


5.00 


± 


0.08 


4.07 


± 


0.08 


5.64 


± 


0.09 


5.53 


± 


0.11 


5.08 


± 


0.20 


6.18 


± 


0.09 








WASP-4 


6.25 


± 


0.12 


3.09 


± 


0.14 


4.87 


± 


0.09 


3.86 


± 


0.06 


5.59 


± 


0.09 


5.53 


± 


0.13 


4.87 


± 


0.06 


6.14 


± 


0.13 


2.14 


± 


0.09 


WASP-5 


6.49 


± 


0.09 


3.32 


± 


0.16 


5.07 


± 


0.12 


4.07 


± 


0.08 


5.77 


± 


0.08 


5.84 


± 


0.09 


5.11 


± 


0.07 


6.37 


± 


0.10 


2.30 


± 


0.14 


WASP-6 


6.19 


± 


0.12 


3.02 


± 


0.10 


4.83 


± 


0.08 


3.80 


± 


0.09 


5.50 


± 


0.10 


5.27 


± 


0.12 


4.67 


± 


0.14 


6.04 


± 


0.12 








WASP-7 














5.26 


± 


0.24 








5.77 


± 


0.06 














6.37 


± 


0.10 








WASP-8 


6.56 


± 


0.08 


3.23 


± 


0.06 


5.14 


± 


0.10 


4.19 


± 


0.09 


5.84 


± 


0.05 


5.66 


± 


0.06 


5.21 


± 


0.17 


6.49 


± 


0.11 


2.37 


± 


0.09 


WASP-15 


6.31 


± 


0.09 


3.17 


± 


0.10 


4.97 


± 


0.07 








5.62 


± 


0.08 


5.30 


± 


0.08 


4.77 


± 


0.10 


6.12 


± 


0.07 


2.29 


± 


0.08 


WASP-16 


6.44 


± 


0.12 


3.14 


± 


0.13 


5.03 


± 


0.10 


4.00 


± 


0.11 


5.73 


± 


0.10 


5.67 


± 


0.16 


5.10 


± 


0.16 


6.41 


± 


0.11 


2.12 


± 


0.05 


WASP-17 


6.27 


± 


0.19 








4.91 


± 


0.06 








5.53 


± 


0.05 


5.19 


± 


0.09 








6.06 


± 


0.10 


2.20 


± 


0.11 


WASP-18 


6.48 


± 


0.12 








5.07 


± 


0.06 


3.94 


± 


0.08 


5.75 


± 


0.10 


5.51 


± 


0.15 








6.24 


± 


0.06 


2.41 


± 


0.10 


WASP-19 


6.51 


± 


0.13 


3.29 


± 


0.09 


5.15 


± 


0.11 


4.16 


± 


0.08 


5.81 


± 


0.08 


5.78 


± 


0.21 


5.12 


± 


0.10 


6.40 


± 


0.08 


2.39 


± 


0.16 



Table 7. Abundance uncertainties for the HARPS solar spec- 
trum. 



Element 


A %s 


A logg 


A l-mic 




+83 K 


+0.11 dcx 


+0.11 kms-1 


[Fe/H] 


0.02 


0.00 


-0.02 


[Ca/H] 


0.05 


-0.03 


-0.02 


[Sc/H] 


0.01 


0.04 


0.00 


[Ti/H] 


0.05 


0.01 


-0.01 


[V/H] 


0.09 


0.00 


0.00 


[Cr/H] 


0.04 


0.01 


-0.01 


[Mn/H] 


0.07 


0.00 


-0.02 


[Co/H] 


0.06 


0.01 


0.00 


[Ni/H] 


0.04 


0.00 


-0.01 


[Y/H] 


0.00 


0.05 


-0.01 


! 8. Abundance uncertainties for Procyon. 


Element 


ATeff 


A logg 


A I'mic 




+83 K 


+0.11 dcx 


+0.11 kms-1 


[Fe/H] 


0.03 


0.01 


-0.01 


[Ca/H] 


0.05 


-0.01 


-0.02 


[Sc/H] 


0.04 


0.04 


0.00 


[Ti/H] 


0.03 


0.02 


-0.01 


[V/H] 


0.05 


0.01 


0.00 


[Cr/H] 


0.02 


0.02 


-0.02 


[Mn/H] 


0.11 


0.06 


0.05 


[Co/H] 


0.05 


0.00 


0.00 


[Ni/H] 


0.05 


0.00 


0.00 


[Y/H] 


0.03 


0.04 


-0.01 



potential. It was also essential that the abundances of Fe i 
and Fe ll agreed. The logg was also determined from the 
Na D and Ca I line at 6439 A. We tested these methods 
on two different solar spectra and a spectrum of Procyon 
before being applied to the WASP stars. The masses and 
radii of the stars, which are needed to obtain planetary pa - 
rameters, were calculated based on the [Torres et al.l (|2010l ) 
calibration. 

We found that the log A{¥e) values determined from the 
HARPS spectra were an average of 0.09 ± 0.05 dex higher 
than those obtained from the lower S/N spectra used in the 
initial analysis of the WASP stars. This is most likely due 
to the higher S/N of the HARPS spectra allowing for the 
wings of the spectral lines to be measured with greater pre- 
cision. We discussed the importance of clearly stating the 
solar log(j4) values. Stellar metallicities are often given rela- 
tive to the Sun, but the solar value used is not always stated, 
which can cause additional discrepancies when comparisons 
are made between different methods. The average uncer- 
tainty in Toff, logp, Vmic and log^(Fe) is 83 K, 0.11 dex, 
0.11 kms^^ and 0.10 dex respectively, showing that there is 
a limit to the accuracy that can be achieved even using high 
S/N spectra. 
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solar atlas, and also from the HARPS Procyon spectrum to 
create a subset of lines for hotter stars. The line list was 
cross referenced with the NIST database in order to reject 
any E-rated lines, i.e. lines with a log g f uncertainty greater 
than ± 50 per cent, however not all lines in the list were also 
present in the database. 

We determined elemental abundances via measuring 
EWs, and a least squares fit of spectral lines was performed 
when the synthetic profile failed to agree with the observed 
spectrum. The continuum was normalised by eye over a 
small wavelength range for each line, however even for high 
S/N spectra it was found that some uncertainty still exists 
in continuum placement. We determined the T^s, logg, and 
Vmic from Fe lines by requiring that there is no trend present 
when abundance is plotted against either EW or excitation 
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